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20.  ABSTRACT  (Continue  on  revetee  eide  II  neceeeery  end  Identity  by  block  number ) 

Energetic  electron  measurements  from  a  particle  detector  on  the  DMSP/F2 
satellite  have  been  used  to  determine  the  corrected  geomagnetic  latitude 
(Ap  )  of  the  equatorward  boundary  of  the  auroral  oval.  The  satellite  was 
launched  into  a  nearly  sun-synchronous,  polar  orbit  centered  on  the  0700-1900 
MLT  meridian.  Due  to  the  wobble  of  the  dipole  and  a  very  slow  precessional 
motion  of  the  orbit  more  than  6000  boundary  crossings  could  be  studied  in  the 
1600-2300  and  0400-1000  MLT  sectors.  The  detectors,  which  have  large  geometric 
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factors,  look  radially  away  from  the  earth  and  detect  precipitating  electrons 
with  energies  between  SOeV  and  20keV.  In  the  evening  sector  the  equatorward 
boundaries  are  precisely  determined  from  the  rise  in  the  total  electron  flux. 
The  morningside  boundary  cannot  always  be  clearly  delineated.  Data  were 
divided  into  one  hour  magnetic  local  time  bins  with  ^400  samples  per  bin.  In 
each  of  the  MLT  bins  was  found  to  be  linearly  correlated  to  the  Kp  index. 

Regression  values  of  were  projected  to  the  magnetic  equatorial  plane 

using  the  Mead-Fairf iela magnetic  field  model.  The  projected  boundaries  are 
not  in  good  agreement  with  the  injection  boundaries  of  Mauk  and  Mcllwain 
(1974) ...  They  are  best  fit  to  boundaries  derived  using  Volland-Stern  type 
electrics,  f ields  with  the  axis  of  symmetry  rotated  in  the  evening  sector  for 
Kp  <  2  and  into  the  afternoon  sector  for  Kp  i  2. 


INTRODUCTION 


The  work  reported  here  is  a  continuation  of  a  larger  effort  to  model 
particle  entry  into  the  extremely  high  latitude  regions  of  the  ionosphere 
(i  80°CGM) .  When  of  sufficient  density  and  energy,  electrons  precipitating 
here  result  in  polar  cap  auroras  (more  generally  referred  to  as  extremely 
high  latitude  auroras) .  A  large  statistical  study  using  DMSP  (Defense 
Meteorological  Satellite  Program)  optical  data  (Gussenhoven,  1979)  showed 
that  the  occurrence  of  auroras  at  these  latitudes  depended  primarily  on  the 
2-component  of  the  interplanetary  magnetic  field  being  positive.  They  also 
occurred  when  magnetic  activity,  as  measured  by  Kp,  was  moderate  (2-3),  while 
electro jet  activity,  as  measured  by  A^,  was  low.  Further,  extremely  high 
latitude  auroras  occur  preferentially  in  local  time  sectors  of  the  cap 
according  to  the  y-component  of  the  interplanetary  magnetic  field,  By.  They 
have  greater  frequency  of  occurrence  in  the  morning  for  negative  B^,  and  in  the 
evening  for  positive  B^.  This  is  in  the  opposite  sense  of  increases  in  the 
polar  cap  electric  field  (Heppner ,  1972). 

Polar  cap  arcs  can  be  intense  and  highly  structured,  but  more  often  they  are 
weak,  and  difficult  to  detect  in  the  DMSP  images.  With  the  launch  of  the 
DMSP-F2  satellite,  in  September,  1977,  which  carried  two  sets  of  electrostatic 
analysers  (J-sensors)  with  large  geometric  factors,  (Hardy  et  al. ,  1979) 
electron  precipitation  in  the  polar  caps  could  be  studied  in  greater  detail. 
Several  characteristics  of  the  electron  precipitation  were  apparent:  1)  Morning 
precipitation  patterns  differ  substantially  from  evening  precipitation  patterns. 
In  the  morning  low  energy  electrons  (<  IkeV)  precipitate  over  much  greater 
latitudinal  ranges  extending  frequently  to  latitudes  >  80°CGlp>>5$pcursion  of 
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electron  precipitation  Into  extremely  high  latitudes  also  occurs  In  the  evening, 
but  much  less  frequently.  That  is,  there  is  not  the  symmetry  In  occurrence 
that  would  be  expected  if  (positive  and  negative  values  occur  at  near 
equal  frequencies)  provided  the  principle  driving  mechanism. 

In  order  to  investigate  the  source  of  particles  precipitating  in  the 
polar  caps  information  on  the  overall  state  of  the  magnetosphere  is  necessary. 

And  in  particular,  it  is  important  to  specify  the  convection  electric  field 
since  it  directly  influences  motion  of  low  energy  particles.  The  inner  edge  of 
plasmasheet  is  a  direct  consequence  of  the  configuration  of  the  magneto- 
spheric  electric  field.  Various  theoretical  models  of  the  convection  electric 
field  exist  in  the  literature,  but  direct  comparison  to  measured  values  is 
nearly  non-existent.  Model  predictions  can  be  more  easily  tested.  Specifically, 
the  DMSP  precipitating  electron  data  give  auroral  equatorward  boundaries 
which  map  into  the  inner  edge  of  the  plasmasheet.  The  position  of  the  boundarv 
is  an  indication  of  magnetospheric  activity,  and  in  particular  of  the  strength 
and  shape  of  the  magnetospheric  convection  electric  field. 

This  report  consists  of  two  completed  studies:  The  first  is  "DMSP-F2 
electron  observations  of  equatorward  auroral  boundaries  and  their  relationship 
to  magnetospheric  electric  fields"  and  is  submitted  to  the  Journal  of  Geophysical 
Research  for  publication.  The  second  is  "Comment  on  'Diurnal  variation  of  the 
auroral  oval  size’  by  C.-I.  Meng"  which  is  currently  in  press  at  the  Journal 
of  Geophysical  Research.  Both  studies  were  a  collaboration  with  D.A.  Hardy 
of  Air  Force  Geophysics  Laboratory,  and  W.  J.  Burke,  at  Boston  College,  and 
are  presented  as  they  were  submitted  for  publication 
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ABSTRACT 


,  ^-Energetic  electron  measurements  from  a  particle  detector  on  the 

DMSP/F2  satellite  have  been  used  to  determine  the  corrected  geomagnetic 

latitude  (ACGM)  of  the  equatorward  boundary  of  the  auroral  oval .  The 

satellite  was  launched  into  a  nearly  sun-synchronous,  polar  orbit  centered 

on  the  0700-1900  MLT  meridian.  Due  to  the  wobble  of  the  dipole  and  a 

very  slow  precessions]  motion  of  the  orbit  more  than  6000  boundary  crossings 

could  be  studied  in  the  1600-2300  and  0400-1000  MLT  sectors.  The  detectors, 

which  have  large  geometric  factors,  look  radially  away  from  the  earth,  and 

detect  precipitating  electrons  with  energies  between  50eV  and  20keV.  In 

the  evening  sector  the  equatorward  boundaries  are  precisely  determined 

from  the  rise  in  the  total  electron  flux.  The  morningside  boundary  cannot 

always  be  clearly  delineated.  Data  were  divided  into  one  hour  magnetic  local 
.  i 

time  bins  with  ^400  samples  per  bin.  In  each  of  the  MLT  bins  .\„™.  was 
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found  to  be  linearly  correlated  to  the  Kp  index.  Regression  values  of 

r  ; 

^CGM  were  projected  to  the  magnetic  equatorial  plane  using  the  Mead- 

r  1 

Fairfield  magnetic  field  model.  The  projected  boundaries  are  not  in  good 
agreement  with  the  injection  boundaries  of  Mauk  and  Mcllwain  (1974). \  They 
are  best  fit  to  boundaries  derived  using  Volland-Stern  type  electric  I  fields 
with  the  axis  of  symmetry  rotated  in  the  evening  sector  for  Kp  <  2  an\l 


into  the  afternoon  sector  for  Kp  2. 
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INTRODUCTION 


The  dynamics  and  sources  of  auroral  electrons  have  long  been  of  concern 
within  the  geophysical  community.  Whether  "inverted-V"  (Frank  and  Acker son. 

1971)  electrons  directly  come  from  the  magnetosheath  (Frank.  1971),  the  magneto- 
spheric  boundary  layer  (Winningham  et  al..  1975),  the  neutral  sheet  (Speiser,  1967), 
or  some  other  location  in  the  magnetosphere  remains  a  controversial  topic.  It 
appears  to  be  widely  accepted,  however,  that  diffuse  auroral  electrons  come 
from  the  central  plasmasheet  (Lui,  et  al..  1977).  Thus,  it  is  expected  that 
the  equatorward  boundary  of  the  auroral  oval  reflects  the  dynamics  of  the 
position  of  the  inner  boundary  of  the  plasmasheet  (Rp) .  In  a  general  sense 
this  appears  to  be  true  in  that  the  equatorward  boundary  of  the  oval  moves 

equatorward  and  Rp  moves  earthward  as  the  level  of  magnetic  activity  increases 

(Vasyliunas .  1968) . 

Studies  of  the  position  of  the  equatorward  boundary  of  the  auroral  oval 

have  been  conducted  using  imagery  from  the  ISIS  2  (Lui,  et  al. ,  1975)  and  DMSP 

(Sheehan  and  Carovillano,  1978)  satellites.  In  both  instances  the  boundary  was 
observed  to  move  equatorward  with  increasing  magnetic  activity.  A  more 
comprehensive  study  by  Kamide  and  Winningham  (1977)  used  data  from  the  soft 
particle  spectrometers  on  ISIS  1  and  2  to  identify  the  position  of  the  boundary 
and  to  show  that  it  shifts  by  ±4°  with  IMF  Bz  +  5nT.  The  statistical  sample 
was  large  enough  to  show  that  for  the  same  value  of  B2,  further  equatorward 
shifts  occur  1)  if  a  substorm  is  in  progress,  2)  during  the  summer,  and  3) 
with  increasing  magnetic  local  time.  A  comparison  of  the  low-latitude  boundary 
of  the  oval  projected  into  the  magnetic  equatorial  plane  with  the  shape  of 
calculated  zero-energy  Alfvdn  layers  (Klvelson,  1976;  Southwood  and  Kaye,  1979) 


met  with  mixed  success.  For  the  remainder  of  this  paper  we  refer  to  the  work 
of  Sheehan  and  Carovillano  (1978)  and  Kamide  and  Winningham  (1977)  as  S-C 


and  K-W,  respectively. 

The  purposes  of  this  report  are  to:  1)  present  comprehensive  measurements 
of  the  corrected  magnetic  latitude  of  the  equatorward  boundary  of  the  auroral 
oval  and  2)  compare  these  results  to  model  predictions  of  Rp,  the  inner 

boundary  of  the  plasmasheet.  The  measurements  were  made  using  an  electron 
spectrometer  on  board  the  DMSP  (Defense  Metereological  Satellite  Program) 
satellite  F2.  In  the  following  section  we  describe  the  instruments,  the 
criteria  for  identifying  A  ,  and  the  methods  used  for  binning  data.  It  is 
shown  that  negative  linear  relationships  exist  between  A  ^  and  Kp  in  both 
evening  and  morning  sectors.  The  discussion  section  is  divided  into  three 
parts:  1)  A  comparison  with  ionospheric  measurements  of  S-C  and  K-W  is 
made.  The  comparison  shows  a  high  level  of  agreement  between  the  three  data 
sets.  2)  The  DMSP  boundaries  are  projected  via  the  Mead-Fairf ield  magnetic 
field  model  to  the  magnetic  equator  and  compared  with  the  injection  boundaries 
of  Mauk  and  Me II wain  (1974).  We  find  the  agreement  to  be  poor.  3)  The 
boundaries  are  then  compared  to  the  semi-empirical  model  of  magnetospher ic 
electric  fields  of  Volland  (1974)  and  Stern  (1975).  The  best  fits  to  DMSP 
observations  are  found  using  shaping  factors  between  2  and  3  and  with  the 
axis  of  symmetry  rotated  away  from  the  dawn-dusk  meridian. 


INSTRUMENTATION  AND  DATA  SELECTION 

DMSP/F2  is  a  three-axis  stabilized  satellite  in  a  near  sun-synchronous, 
circular  orbit  at  an  altitude  of  840  km.  The  orbital  period  is  101  minutes; 
the  nominal  inclination  is  98.75°.  At  launch  the  orbit  was  centered  near  the 
0700-1900  meridian  but  is  subject  to  a  very  slow  precession  toward  later  local 
times.  Due  to  the  offset  between  the  earth  spin  axis  and  magnetic  axis,  the 
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orbit  is  subject  to  significant  diurnal  and  seasonal  variations  in  the  magnetic 
local  time-magnetic  latitude  frame  of  reference.  Thus,  A„_.  is  sampled  in  a 
wider  range  of  MLT  than  might  be  assumed  from  the  restricted  geographic  local 
time  locations  of  the  orbit.  It  should  also  be  noted  that  DMSP  is  an  operational 
Air  Force  satellite.  Except  during  periods  of  down-link  transmissions,  data 
are  almost  always  being  recorded. 

The  particle  detector  on  DMSP/F2  consists  of  two  curved  plate  electro¬ 
static  analyzers  that  measure  the  fluxes  of  electrons  in  16  energy  channels 
between  50eV  and  20keV  once  per  second.  The  apertures  of  the  analyzers  always 
face  local  vertical  such  that  at  auroral  and  polar  cap  latitudes  they  detect 
precipitating  rather  than  back  scattered  and/or  trapped  electrons.  One 
analyzer  covers  the  energy  range  from  50eV  to  IkeV  with  a  geometric  factor  of 
4  x  10-Ucm2  ster  and  a  AE/E  of  10  percent.  The  other  analyzer  covers  the 
energy  range  from  IkeV  to  20keV  with  a  geometric  factor  of  10-3cm2  ster  and  a 
AE/E  of  12  percent.  The  large  geometric  factors  insure  that  the  flux  level 
for  the  electrons  in  the  diffuse  aurora  is  well  above  the  detector's  sensitivity. 
A  detailed  description  of  the  detector  is  given  by  Hardy  et  al.  (1979). 

An  example  of  DMSP/F2  electron  data,  taken  from  a  south  polar  pass  on 
2  July,  1977  is  given  in  Figure  1.  Data  are  plotted  as  JTOT,  the  directional 
integral  flux  (cm2ster  sec)-1  in  the  bottom  panel;  JETOT,  the  directional 
energy  flux  (keV/cm2ster  sec)  in  the  middle  panel,  and  EAVE,  the  average 
energy  in  keV  in  the  top  panel.  The  scale  for  EAVE  is  linear.  These  quantities 
are  plotted  as  functions  of  universal  time  in  seconds  of  the  day;  the  geographic 
latitude  and  longitude  of  the  subsatellite  position;  the  corrected  geomagnetic 
latitude  and  longitude  of  the  satellite  projected  to  an  altitude  of  100  km; 
and  the  magnetic  local  time.  Values  of  the  equatorial  boundaries  (A  )  of 

ton 

the  auroral  oval  are  assigned  to  the  corrected  geomagnetic  latitudes  >45° 


at  which  JTOT  rises  noticeably  above  background.  This  is  essentially  the  same 
method  used  by  Kamide  and  Winninghara  (1977).  In  the  example  given  in  Figure  1, 
J\CGM  is  at  -66.6°  on  the  evening  side  (^20MLT)  and  at  -65.1°  on  the  morning 
side  0^07MLT). 

Several  characteristics  of  the  electron  precipitation  that  pertain  to 

the  choice  of  auroral  boundaries  are  illustrated  in  Figure  1.  First,  we 

have  found  the  condition  JTOT>107  (cm2ster  sec)-1,  indicated  in  the  figure 

by  the  horizontal  line,  provides  a  generally  useful  criterion  for  selecting 

both  the  poleward  and  equatorward  boundaries  of  the  oval.  The  poleward 

boundaries  selected  in  the  case  given  in  Figure  1  are  marked  for  purposes 

of  illustration.  They  and  their  attendant  selection  problems  will  not  be 

discussed  here.  Second,  the  evening  equatorward  boundary  is  sharper  than 

its  morningside  counterpart.  On  the  morning  side  JTOT  rose  from  background 

to  107  (cm2ster  sec)-1  over  1.4°  latitude,  whereas  the  evening  side  rise  was 

nearly  instantaneous.  Generally,  flux  gradients  are  found  near  both  equatorvard 

boundaries.  The  gradient  never  extends  more  than  a  few  degrees  latitude  in  the 

evening  sector,  but  can  cover  as  much  as  10°  on  the  morning  side  where  it  is 

always  perceptibly  present.  Third,  equatorward  of  the  evening  and  morning 

values  of  h  there  are  slight  rises  in  JTOT  with  relatively  large  values  of 
COM 

JETOT  and  EAVE.  These  increases  are  due  to  radiation  belt  particles  that 
penetrate  the  detector  casing  and  directly  stimulate  the  channeltrons ;  they 
must  be  differentiated  from  the  auroral  electrons  when  determining  boundaries. 
Due  to  the  geometrical  configurations  of  the  two  detectors,  radiation  belt 
contamination  is  largely  limited  to  the  energy  channels  of  the  1  to  20  keV 
detector. 

The  boundary  of  the  oval,  ^CGM»  cannot  always  be  determined  unambiguously. 
On  both  morning  and  evening  sides  of  the  oval,  the  onset  of  electron  precipi¬ 
tation  can  be  obscured  in  various  ways.  In  the  evening  the  principle  problem 
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is  contamination  of  the  low  energy  channels  of  the  detector  by  photoelectrons. 

As  an  example  in  Figure  2  we  have  plotted  the  count  rates  in  representative 
channels  of  the  detector  as  a  function  of  increasing  geographic  latitude. 

Arrows  show  where  there  is  a  significant  increase  in  counts  due  to  the  onset  of 
precipitating  auroral  electrons.  The  increase  in  counts  is  observed  first  in 
the  lowest  energy  channel  and  then,  progressively,  in  higher  energy  channels. 

In  the  50eV  channel  significant  counts  are  observed  equatorward  of  the  arrow. 
These  counts  arise  from  photoelectrons  from  the  conjugate  ionosphere.  The  effect 
of  photoelectrons  can  be  observed  at  energies  as  high  as  several  hundred  electron 
volts.  Since  the  photoelectrons  are  observed  in  the  same  channels  as  the  onset 
of  precipitation,  the  ambiguity  in  the  boundary  is  introduced  by  having  to 
separate  the  two  effects.  Also  shown  in  Figure  2  are  increased  counts  in 
the  5.5keV  channel  equatorward  of  auroral  electron  precipitation.  This  is  an 
example  of  radiation  belt  contamination.  As  discussed  above  the  effect  is 
seen  only  in  the  channels  of  the  high  energy  detector.  In  the  evening 
sector,  auroral  precipitation  is  first  detected  in  the  lowest  energy  channels. 

At  times  when  the  radiation  belt  and  the  auroral  oval  overlap  the  combination 
of  high  energy  channel  contamination  and  increased  low  energy  count  rates 
produces  an  increase  in  JTOT  and  a  decrease  in  EAVE.  Thus,  the  auroral 
boundary  determination  is  unaffected. 

Uncertainties  in  identifying  the  position  of  the  boundary  in  the  morning 
sector  are  more  severe.  There  are  three  sources  of  uncertainty:  1)  overlap 
between  regions  of  energetic  electron  precipitation  and  radiation  belt 
contamination  of  corresponding  energy  channels;  2)  a  much  more  gradual 
latitudinal  onset  of  precipitation  than  is  found  in  the  evening,  and  3)  the 
existence  of  energetic  plasma  close  to  but  detached  from  the  boundary.  In 


< 


3 


Figure  3  we  have  plotted  for  the  morning  auroral  zone  the  counts  as  in  Figure  2 
The  figure  shows  that  for  this  case  the  onset  of  electron  precipitation 
in  the  morning  appears  first  in  the  high  energy  channels  and  then  progressively 
in  lower  energy  channels.  For  this  case  the  region  in  which  penetrating 
particles  are  producing  counts  are  equatorward  of  the  onset  of  auroral 
precipitation.  In  active  periods  they  can  overlap.  Since  the  count  levels 
are  comparable,  for  the  cases  of  overlap  it  is  difficult  to  separate  the  two 
effects  and  our  ability  to  determine  the  boundary  is  affected. 

The  second  effect  causing  ambiguity  in  identifying  the  morningside  boundar 

is  illustrated  in  Figure  4  which  shows  a  complete  south  polar  pass  from 

28  February,  1978.  On  the  morning  side  a  gradient  in  JTOT  continues  into  the 

region  of  radiation  belt  contamination.  JETOT  and  EAVE  fall  sharply  at 

-64.7°CGM  but  JTOT  remained  too  high  to  be  explained  as  contamination  until 

-58.5°.  At  -58.5°  EAVE  rose  with  significant  scatter  indicating  the  presence 

of  dominant  radiation  belt  contamination  and  obscuring  the  possible  presence 

of  keV  electron  precipitation.  In  cases  such  as  this  the  boundary  was  chosen 

to  be  at  the  lowest  magnetic  latitude  at  which  precipitating  electrons  could 

be  identified;  in  this  case  A__.,  =  -58.3°. 

•  t/GM 

The  third  effect  which  can  be  found  in  both  morning  and  evening  regions 
is  shown  in  Figure  5.  This  example  shows  a  north  pole  pass  on  January  21,  1°77 
during  a  time  of  magnetospheric  quiet.  On  the  evening  side  there  are  small 
regions  of  energetic  electrons  near  to  but  detached  from  the  more  poleward 
precipitation  region.  These  can  occur  during  active  times,  as  well,  and  are 
Ignored  in  selecting  the  boundaries. 

The  above  effects  produced  greater  than  1  degree  uncertainty  in  the 
equatorial  boundary  only  for  a  small  percentage  of  the  total  cases.  On  the 


evening  side  approximately  90  percent  of  boundaries  were  determined  with 
no  ambiguity.  On  the  morning  side  the  figure  was  approximately  70  percent. 

More  than  6000  auroral  oval  boundary  crossings  were  determined  using  data 
from  September,  1977  to  December,  1978.  Over  two-thirds  of  the  crossings 
occurred  during  north  pole  Fall  and  Winter  months.  It  was  necessary  to  use 
Spring  and  Summer  data  to  obtain  a  good  statistical  sample  for  high  Kp  values. 
The  measurements  are  binned  into  one  hour  units  of  MLT  and  to  all  steps  in 
Kp  (0,  0+,  1  ...).  Cases  were  selected  to  insure  that  20  to  25  samples 

were  included  in  each  Kp  <  6  element.  In  some  MLT  sectors  a  statistically 
significant  number  of  samples  were  available  for  Kp  >  6  elements.  Due  to 
the  orbital  characteristics  of  DMSP/F2  these  criteria  were  largely  met  in 
the  17-23  and  04-10  MLT  sectors  of  the  northern  hemisphere  and  the  17-21 
and  04-08  sectors  of  the  southern  hemisphere.  The  major  exceptions  were  the 
18-19  MLT  sector  of  the  northern  hemisphere  in  which  only  103  samples  were 
found,  and  the  15-16  MLT  sector  in  the  southern  hemisphere  having  107  samples. 
The  small  sample  size  resulted  from  the  fact  that  the  DMSP/F2  orbit  drifted 
slowly  out  of  (into)  these  regions  over  the  time  interval  used. 

OBSERVATIONS 

Examples  from  the  evening  (19-20)  and  morning  (07-08)  MLT  sectors  in 
which  coverage  is  available  from  both  the  northern  (top  traces)  and  southern 
(bottom  traces)  hemispheres  are  given  in  Figures  6  and  7,  The  average  value 
and  the  standard  deviation  of  the  boundaries  within  each  Kp  bin  is  shown. 

Solid  lines  represent  linear  regressions  perfo’"  ed  with  all  the  points  in 
the  MLT  sector  versus  Kp  (not  with  the  mean  values).  The  equation  for  the 
line  is  given  in  the  upper  right  hand  corner  of  each  plot.  The  correlation 
coefficient  (cc)  is  also  given. 
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From  an  examination  of  data  contained  in  these  figures  several  points 
can  be  made:  1)  Kp  orders  the  equatorvard  auroral  boundaries  linearly 
throughout  the  entire  range.  The  mean  values  lie  on  or  very  near  the  line. 

The  fit  is  exceptional  for  the  evening  values,  and  while  there  is  more  varia¬ 
tion  in  the  morning,  no  trend  other  than  linear  is  apparent.  The  correlation 
coefficients  across  the  whole  Kp  range  are  very  high:  <  -.81  for  the  four 
cases  shown.  2)  Standard  deviations  about  mean  values  tend  to  increase  with 
increasing  Kp.  Also,  standard  deviations  in  the  morning  sector  are  greater 
than  those  found  in  the  evening  sector.  This  reflects  our  general  observation 
that  auroral  electron  precipitation  patterns  are  far  more  variable  in  the 

morning  than  evening  sector.  3)  In  MLT  regions  of  overlap  A__,  has  lower 

CGM 

values  in  the  southern  hemisphere.  The  difference  is  greater  in  the  morning 
for  low  Kp  values.  The  slopes  in  the  straight  line  fit  have  larger  negative 
values  in  the  north  pole,  and  the  difference  between  the  boundaries  approaches 
zero  for  large  Kp.  The  difference  may  be  due  to  north-south  assymmetries 
in  mirror  points  of  plasmasheet  electrons.  It  may  also  reflect  a  summer- 
winter  assymmetry  due  to  the  tilt  of  the  dipole  noted  by  Kamide  and  Winningham 
(1977).  Recall  that  most  of  the  data  points  used  in  the  present  study  were 
acquired  within  three  months  of  the  December  1977  solstice. 

The  entire  data  set  is  summarized  for  the  morning  and  evening  sectors  in 
Tables  1  and  2,  respectively.  Here  we  have  listed  values  of  Aq  and  a 
calculated  from  the  regression  formula:  A M  -  ±(A  +  aKp) ,  where  the  +(-) 

sign  refers  to  the  north  (south)  pole.  Also  listed  are  the  number  of  cases, 

N,  and  the  correlation  coefficient,  cc.  Even  including  the  MLT  sectors  for 
which  the  statistical  sample  was  small  the  correlation  coefficients  lie  in 
the  range  -.58  >  cc  >  -.90.  For  the  full  samples  correlation  coefficients 
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were  less  Chan  -.7  in  the  morning  and  less  than  -.8  in  the  evening.  This 
suggests  that  Kp  is  better  for  organizing  than  IMF  Bz-  Correlation 

coefficients  reported  by  K-W  were  in  the  -0.4  >  cc  >  -0.7  range.  Evening 
sector  values  of  A^,^,  for  a  given  value  of  Kp  decreased  with  increasing 
MLT.  According  to  K-W  this  decrease  continues  into  the  post-midnight  sector. 
In  the  dawn  sector  increases  with  increasing  MLT.  However,  at  Kp  ■  6 
dawn  and  post-dawn  sector  values  of  approach  a  constant  value  of  '-57°. 

Values  of  a  in  the  evening  sector  fall  in  the  restricted  range  -1.6  i  0.4, 
and  -1.75  i  0.4  in  the  morning  sector.  The  largest  negative  values  of  a 
were  at  07-08  MLT  and  18-19  MLT,  indicating  that  the  greatest  variability  in 
the  boundaries  occurs  near  dawn  and  dusk. 

DISCUSSION 

A.  COMPARISON  WITH  MEASURED  IONOSPHERIC  BOUNDARIES 

Direct  quantitative  comparisons  of  results  presented  in  the  previous 
section  with  the  work  of  S-C  and  K-W  are  not  possible  for  many  reasons. 
Correlated  quantities  were  not  the  same  and  were  sampled  largely  in  different 
MLT  sectors.  S-C  used  optical  rather  than  particle  signatures  to  identify 
the  corrected  magnetic  latitude  of  the  equatorward  boundary  as  a  function  of 


Kp.  K-W  studied  the  invariant  latitude  of  the  boundary  as  a  function  of  IMF 
B^ .  Although  Bz  influences  geomagnetic  activity,  its  relationship  with  Kp  is 
not  a  simple  one.  However,  some  semi-qualitative  comparisons  of  the  various 
results  can  be  made  if  we  assume  that  times  when  Bz  E  +(-)5  correspond  to 
periods  of  low  (high)  geomagnetic  activity.  With  these  simplifying  assumptions 


X  refer  to  DMSP/F2  northern  and  southern  hemisphere  observations,  respectively. 

The  high  (low)  latitude  points  were  calculated  for  cases  of  Kp  *  0  (5)  using 

regression  coefficients  of  Tables  1  and  2.  The  A  symbols  refer  to  boundary 

values  calculated  using  regression  coefficients  given  in  Table  1  of  K-W  with 

■  +5y.  The  s  symbols  refer  to  measurements  at  MLT  ■  01,  20,  and  22  for 

Kp  *  0  and  5  as  taken  from  Figure  2  and  8  of  S-C. 

The  auroral  boundary  derived  from  the  three  data  sets  for  periods  of 

magnetic  quiet  shows  a  remarkable  consistency.  In  MLT  sectors  of  common 

measurements  the  quiet  time  results  are  within  1°  of  each  other.  The  boundary 

appears  to  be  continuous  across  the  entire  night  side.  The  one  anomalous  point 

in  the  03-04  MLT  sector  was  based  on  a  data  set  of  8  points  and  is  virtually 

without  statistical  signficance.  It  is  interesting  to  note  that  Meng  et  al. 

(1977)  have  shown  that  the  poleward  boundary  of  the  quiet  auroral  oval  can  be 

represented  by  an  offset  circle.  During  magnetically  quiet  times  the  oval 

tends  to  be  of  narrow  latitudinal  width.  Perhaps  not  surprisingly  A^.  can  also 

be  well  fit  to  an  offset  circle.  (The  symbol  A  is  used  in  comparing  DMSP  with 

r. 

equatorward  boundaries  of  S-C  and  K-W  rather  than  the  symbol  A  which  is 

CUM 

reserved  for  DMSP/F2  data  alone.)  In  Figure  8  the  quiet  time  boundary  is  fit 
to  a  circle  of  radius  21°  that  is  centered  at  88.0°  and  02.7  MLT  (dashed  line). 
DMSP/F2  values  fall  within  the  circle  at  18-20  MLT  and  06-08  MLT,  deviations 
which  will  later  be  shown  to  have  significance. 

At  times  of  high  geomagnetic  activity  agreement  between  the  data  sets  is 
not  as  good.  The  boundary  for  Kp  ■  5  derived  from  DMSP/F2  measurements  is 
•mooth  in  the  evening  and  morning  sectors.  The  high  activity  boundaries  from 
S-C  and  K-W  are  mutually  consistent  but  are  displaced  poleward  of  F2  results 
by  2°  in  the  pre-midnight  region.  No  comparison  is  possible  for  the  post-midnight 
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region.  Variation  from  S-C  data  can  be  accounted  for  as  a  difference  between 
optical  and  particle  boundaries.  On  the  other  hand,  in  their  analysis  K-W 


noted  that  for  the  same  value  of  B  the  value  of  A  could  be  shifted  equator- 

ward  by  several  degrees  if  a  substorm  was  in  progress.  During  times  of 

Kp  ■  5  it  is  expected  that  substorms  would  almost  always  be  in  progress. 

Examining  the  distribution  of  points  in  Figure  3a  of  K-W,  we  see  that  in  the 

20-21  MLT  sector  for  Bz  ■  -5nT,  Ap  was  observed  at  61°  and  59°  where  the 

regression  value  is  63°.  In  general,  most  of  K-W's  data  points  taken  during 

times  when  substorms  were  in  progress  lie  well  below  the  regression  line. 

Evidently,  the  value  of  A  reflects  the  dynamic  state  of  the  magnetosphere 

which  depends  on  physical  parameters  other  than  just  the  hourly  average  of 

B  . 
z 

B.  COMPARISON  WITH  MAGNETOSPHERIC  BOUNDARIES 

The  notion  of  an  injection  boundary  was  first  proposed  by  Mauk  and 
Mcllwain  (1974)  to  systematize  magnetic  local  times  (4>)  at  which  the  ATS-5 
satellite  at  6.7Rp  encountered  sharp  increases  in  low  energy  (keV)  plasma 
fluxes  as  a  function  of  Kp.  A  least  squares  regression  analysis  gave  the 
relationship 

♦  <L  -  6.7 )  -  25.5  -  1.5  Kp  (Kp  £  5+)  .  (1) 

To  extend  these  points  and  to  give  a  generalized  plasma  boundary  Mauk  and 
Mcllwain  used  Carpenter  and  Park's  (1973)  relationship  between  the  position 
of  the  plasmapause  at  $  *  0600  (L^)  and  Kpj  Lp  ■  5.7  -  .47  Kp.  The  Carpenter 
and  Park  relationship  was  derived  using  the  highest  values  of  Kp  in  the 
previous  24  hours  rather  than  its  current  value.  Mauk  and  Mcllwain  assumed 
that  the  shape  of  the  boundary  was  of  the  form 
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(2) 


L(*’  V  "  (TTC)  ‘ 

By  forcing  the  boundary  to  pass  through  L  *  6.7  at  local  times  specified  in 
(1)  it  can  be  shown  that  the  constants  S  and  C  have  values  of  21.4  and  -7.3, 
respectively. 

The  term  injection  boundary  is  not  without  ambiguity  (Kivelson  et  al. 

1978)  .  It  may  either  signify  a  dynamic  plasma  boundary  which  is  suddenly  pushed 
earthward  to  the  position  of  the  satellite  or  refer  to  a  quasi-static  boundary 
that  is  crossed  by  the  satellite.  In  the  dynamic  case  the  electron  flux 
shows  little  or  no  energy-time  dispersion.  In  the  quasi-static  case  electrons 
are  first  detected  at  low  energies,  then  at  higher  energies  as  the  satellite 
penetrates  further  into  the  boundary.  Kivelson  et  al.  (1978)  have  suggested 
that  this  energy-time  dispersion  is  characteristic  of  the  satellite  crossing 
the  Alfvgn  boundaries  of  electrons  with  higher  and  higher  energies.  As 
discussed  below,  Alv£n  boundaries  mark  boundaries  between  open  and  closed 
magnetospheric  trajectories  for  particles  of  given  charge  and  magnetic  moment. 
The  zero  energy  Alfvgn  layer  is  thus  the  boundary  between  open  and  closed 
equipotential  surfaces  in  the  magnetosphere.  For  the  remainder  of  this 
paper,  we  treat  quasi-static  injection  boundaries  and  Alfven  layers  as 
equivalent  concepts. 

The  energy-latitude  dispersion  shown  in  Figure  2  to  be  characteristic 
of  the  equatorward  edge  of  auroral  precipitation  in  the  evening  sector  is 
precisely  what  would  be  expected  from  mapping  a  quasi-static  injection  boundary 
into  the  ionosphere,  and  gives  credance  to  our  association  of  the  DMSP/F2 
boundaries  with  the  inner  edge  of  the  established  plasma  sheet.  A  quasi¬ 
static  injection  boundary  Rp  can  be  determined  from  by  mapping  these 
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latitudes  Into  the  equatorial  plane  where  comparison  can  then  be  made  to 
the  Mauk-Mcllvaln  Injection  boundary. 

Before  making  this  comparison,  several  points  should  be  kept  in  mind. 
First,  the  mapping  Is  extremely  sensitive  to  magnetic  field  models,  and  small 
variations  in  A appear  as  large  variations  in  L-shell.  Figure  9  illustrates 
both  of  these  points.  Here  we  have  mapped  the  best  circle  fit  to  the  quiet 
auroral  oval  shown  in  Figure  8  onto  the  equatorial  plane  using  a)  a  dipolar 
projection  (dashed  line)  and  b)  a  projection  using  the  Mead-Fairf ield  magnetic 
field  model  for  Kp  <  2  (solid  line).  (See  Fairfield  and  Mead.  1975).  The 
difference  between  the  two  can  be  as  large  as  2R^  in  the  evening  sector.  Also 
shown  in  Figure  9  is  the  best  smooth  fit  to  the  data  points,  A  mapped  using 
the  Mead-Fairf ield  magnetic  field  (dash-dot  line).  The  north  and  south  pole 
values  are  averaged  for  overlapping  MLT  sectors.  The  small  deviation  of  ’ 

CO  .i 

from  the  circular  fit  in  the  18-20  MLT  sector,  which  was  pointed  out  previously 

results  in  a  significant  boundary  change  in  the  equatorial  plane. 

Because  the  representation  of  the  Mauk-Mcl lwain  boundary  given  in  (2)  is 

valid  only  at  local  times  prior  to  0200  comparisons  with  DMSP/F2  boundaries  are 

limited  to  the  evening  sector.  We  can,  however,  directly  compare  the  DMSP/F2 

boundary  with  Carpenter  and  Park’s  (1973)  plasmapause  measurements  at  dawn. 

Figure  10  shows  this  comparison.  The  DMSP  values  are  calculated  from  averages 

of  A__,.  in  the  05-06  and  06-07  MLT  sectors.  The  inner  edge  of  the  plasmasheet 
CCM 

as  determined  by  DMSP/F2  auroral  boundaries  lies  beyond  the  plasmapause, 
asymptotically  approaching  a  constant  separation  of  ^0.3Rp  for  large  Kp. 

This  separation  is  within  that  of  the  standard  deviation  in  A  ^  .  The 
deviation  from  linearity  for  small  Kp  cannot  be  accounted  for  in  this  way 
and  indicates  that  the  assumption:  L  *  L^  in  (2)  for  a  given  local  time,  is 
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a  poor  one  for  low 


Figure  11  compares  the  Mauk-Mcllwain  injection  boundaries  (dashed  lines) 


for  Kp  ■  0,  2  and  5,  to  the  inner  edge  of  the  plasmasheet  as  determined 
from  values  for  the  same  Kp  levels  (solid  lines).  For  Kp  2  and  5  the 

Mauk-Mcllwain  injection  boundaries  lie  within  the  plasmasheet.  This  is 
consistent  with  photometric  auroral  observations  at  the  foot  of  the  ATS5 
field  line  by  Eather  et  al.  (1976).  They  found  that  substorm  related  injections 
occurred  within  the  existing  plasmasheet  in  the  pre-midnight  sector.  These 
injections  were  furthest  from  Rp  at  dusk  and  approached  Rp  near  midnight. 
Mauk-Mcllwain  statistics  include  all  injection  boundaries  encountered,  a  mix 
of  dynamic  and  quasi-static  boundaries.  The  boundary  that  results  (R.,^  ) 
should  then  lie  within  the  plasmasheet,  or  at  R^  >  Rp-  This  is  the 
relationship  the  Mauk-Mcllwain  boundary  has  to  the  DMSP/F2-  determined  boundary 
for  Kp  «  2-5;  with  R^^  ■+  Rp  near  midnight. 

For  Kp  *  0,  the  two  boundaries  differ  greatly,  particularly  in  the  evening 
sector.  The  local  time  at  which  ATS-5  statistically  crosses  the  quiet-time 
injection  boundary  is  0130,  which  point  compares  very  favorably  with  the  inner 
edge  of  the  plasmasheet  as  measured  by  DMSP/F2.  Using  higher  values  for  the 
inner  edge  of  the  plasmasheet  at  dawn  than  those  given  by  Carpenter  and  Park 
will  expand  the  Mauk-Mcllwain  boundary,  but  the  essentially  different  shape  of 
the  DMSP/F2  boundary  indicates  that  for  low  Kp  values  the  Mauk-Mcllwain  boundary 
is  valid  over  a  very  restricted  local  time  range.  This  is  not  surprising  since 
geosynchronous  orbits  rarely  enter  the  plasma  sheet  during  very  quiet  times, 
and  when  they  do,  it  is  post-midnight,  far  from  the  dusk  boundary  in  question. 


C.  A  QUASI -EMPIRICAL  REPRESENTATION  OF  THE  OBSERVED  BOUNDARIES 

Charged  particles  in  the  magnetosphere  are  subject  to  gradient-curvature 
and  to  E  x  B  drifts.  The  gradient-curvature  drift  depends  on  the  charge  and 
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energy  of  the  particle  while  the  electric  field  drift  does  not.  The  electric 
field  is  a  superposition  of  the  general  dawn  to  dusk  convective  field,  the 
corotation  field  and  fields  due  to  charge  separations  within  the  system.  The 
latter  type  of  field  is  primarily  due  to  shielding  charges  near  the  plasma- 
pause  (Jaggi  and  Wolf,  1973)  and  to  polarization  charges  at  conductivity 
gradients  near  the  terminator  in  the  ionosphere  (Wolf ,  1970).  The  shape  and 
position  of  magnetospheric  boundaries  for  low  energy  plasmasheet  electrons 
under  the  assumption  of  a  uniform  cross-tail  electric  field  has  been  discussed 
by  Kivelson  (1976).  The  electrons  drift  earthward  from  the  tail  along  equi- 
potential  surfaces,  then  eventually  exit  the  magnetosphere  across  the  dayside 
magnetopause.  The  zero  energy  Alfven  boundary  delineates  regions  of  allowed 
and  forbidden  access.  Its  characteristic  teardrop  shape  is  that  of  the  last 
closed  equipotential .  Based  on  electric  field  measurements  at  ionospheric 
altitudes  (Heppner  1972),  Volland  (1973)  and  Stern  (1975)  suggested  that  the 
convective  field  is  not  uniformly  dawn  to  dusk.  Ejiri  et  al .  (1978)  derived 
an  analytic  expression  for  the  zero  energy  Alfven  boundary  under  nonuniform 
electric  field  situations  that  uses  a  shaping  factor  y.  Recently  Rave  (1980) 
suggested  that  the  value  of  y  reflects  the  degree  to  which  space  charges 
shield  the  convective  electric  field  from  the  inner  magnetosphere.  The  case 
Y  *  1  corresponds  to  situations  of  uniform  dawn  to  dusk  fields  in  which  there 
is  no  shielding.  The  major  effect  of  the  near-terminator  conductivity  gradient 
is  to  rotate  the  dusk  bulge  of  the  Alfv4n  boundary  into  the  evening  sector 
(Wolf,  1970). 

In  the  equatorial  plane  of  the  magnetosphere  the  potential  from 
which  the  total  electric  field  is  derived  is  assumed  to  be  of  the  form: 


nBnV  Y 
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where  fi  is  the  angular  spin  velocity  of  the  earth,  Bq  the  magnetic  field 
strength  at  the  magnetic  equator  on  the  earth's  surface,  R  is  the  distance 
from  the  center  of  the  earth  and  ij>  is  the  local  time  measured  counter¬ 
clockwise  from  midnight.  The  phase  angle,  <j>o>  represents  the  rotation  of 
the  axis  of  symmetry  away  from  the  dawn-dusk  meridian.  The  boundary  between 
open  and  closed  equipotentials  is  calculated  by  assuming  that  at  a  distance 
Rg  and  a  local  time  $  -  4>0  =  3n/2,  the  convection  and  corotational  electric 

fields  cancel  each  other.  The  condition  that  34/3R  =  0  at  this  point  allows 

fiB  Rr3 
o  e. 

a  determination  of  the  constant  A  = - (~Y+T)  '  following  Ejiri  et  al. 

YRS 

(1978)  it  is  easily  shown  that  the  equation  for  the  boundary  is: 
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sin($  -  )  +  (y  +  1)  r - y  =  0 

o  Rg 


Southwood  and  Kaye  (1979)  have  shown  that  an  excellent  approximation  to  the 


solution  of  (A)  is 


R  fl  +  ^ 
S\  1 


r  • 


90  -  (4  -  <p 


where  s  55  cos 


For  y  =  1  the  solution  is  exact.  The 


boundary  has  a  tear-drop  shape,  becoming  less  eccentric  for  higher  values  of  . 

From  the  shape  of  the  boundary  illustrated  in  Figure  9  it  is  obvious 
that  the  axis  of  symmetry  for  the  equatorial  projection  of  the  DMSP/F2  measured 
boundary  is  not  the  dawn-dusk  meridian;  i.e.  ^  4  0.  The  degree  of  rotation 
from  the  dawn-dusk  meridian  as  well  as  the  shapes  of  the  boundary  in  the 
0  <  Kp  <  5  range  are  given  in  Figure  12.  For  Kp  =  0  and  1  the  rotation  is 
counterclockwise  by  23°  and  8°,  respectively.  For  Kp  *  2,  3  and  5  the  axis 
of  symmetry  is  rotated  clockwise  (into  the  afternoon  local  time  sector)  by 
approximately  30°,  A5°  and  45°  respectively.  The  exact  degree  of  clockwise 
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rotation  cannot  be  determined  at  this  time  because  of  the  absence  of  DMSP/F2 
data  in  the  afternoon  and  postmidnight  sectors.  The  rotation  of  the  boundary's 
axis  of  symmetry  into  the  evening  sector  at  times  of  low  Kp  is  due  to 
ionospheric  charge  pile  up  near  the  day-night  terminator  (Wolf ,  1970). 

Increasing  particle  precipitation  at  times  of  high  Kp  weakens  the  ionospheric 
conductivity  gradient  near  the  terminator.  It  is  thus  expected  that  for 
higher  Kp  the  axis  would  rotate  back  toward  the  dawn-dusk  meridian.  The  physical 
basis  for  a  strong  rotation  of  the  symmetry  axis  into  the  afternoon  sector  for 
Kp  1  2  is  not  yet  understood. 

Figure  12  also  shows  results  of  fitting  the  zero  energy  Alfven  boundary 
for  y  =  2  and  3  using  empirical  values  of  £  .  The  method  used  for  deriving 
the  boundary  values  is  as  follows:  1)  Data  points  marked  by  the  symbol  9 
are  obtained  for  a  given  Kp  using  regression  coefficients  from  Tables  1  and 
2.  Where  north  and  south  data  sets  overlap  in  local  time,  average  boundary 
values  are  used.  The  Mead-Fairf ield  model  is  then  used  to  map  points  into  the 
equatorial  plane.  2)  Because  the  boundary  represented  by  (4)  gives  a  shape 
but  not  a  distance  for  a  given  value  of  y,  the  data  points  closest  to  the 
boundary  minimum  (i.e.  opposite  the  dusk  side  stagnation  point)  are  used  to 
'anchor'  the  boundary.  This  procedure  is  justified  because  the  radial  distance 
of  the  boundary  varies  slowly  in  the  morning  local  time  sector  and  the  boundary 
shape  is  nearly  the  same  for  all  values  of  y.  Boundaries  derived  for  the 
y  «  1  case  are  not  shown  in  Figure  12  because  their  stagnation  points  lie 
well  outside  the  y  m  2  values.  Dashed  lines  representing  the  Mauk-Mcllwain 
boundaries  are  provided  for  reference.  The  data  points  in  Figure  12  lie  close 
to  the  y  *  2  and  y  m  3  curves  with  a  tendency  toward  slightly  lower  values  of 
y  for  higher  values  of  Kp.  Values  of  y  between  2  and  3  are  in  agreement  with 
a  variety  of  equatorial  observations  reported  by  Ejiri  et  al.  (1978),  Southwood 
and  Kaye  (1979),  and  Hughes  et  al.  (1979).  Kaye  (1980)  suggested  that  the 
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tendency  toward  lower  values  of  y  with  increasing  Kp  reflects  a  more  dynamic 
mangetosphere  in  which  shielding  layers  do  not  have  sufficient  time  to  become 
fully  developed.  This  hypothesis  awaits  DMSP  data  acquisition  in  the  afternoon 
sector  for  its  verification. 

Data  presented  in  Figure  12  can  also  be  used  to  provide  a  useful,  empirical 
expression  for  the  electrostatic  potential.  Equation  (3)  can  be  written: 


4(L,<p) 


sinU 


-  *o> 


(6) 


where  L  =  R/R^  and  Lg  =  Rg/R^.  It  is  useful  to  define  a  parameter 
Y+l 

£  =  1/Lg  Y .  Using  stagnation  distances  derived  from  Figure  12  as  a  function 
of  L  a  linear  relationship 


£  =  (3.9  +  6. ORp)  •  10"s  (cc  =  0.97),  (7) 

was  found  for  the  case  y  =  2.  In  the  case  y  =  3,  the  quadratic  relationship 


£  =  (8.0  +  4.4  Kp):  •  10“-  (cc  =  0.99)  , 


(S) 


was  found.  Thus  for  the  case  y  =  2  the  total  potential  is 


*(L,$)  =  F.B  R  2 
o  E 


(3.9  -  6 . OKp)  •  10"5  L2sin(f  -  $o)  -  j- 


(9) 


It  should  be  recalled  that  <f>o  is  a  function  of  Kp  that  is  not  describable 
in  a  simple  analytic  form. 


SUMMARY 


By  associating  the  equatorward  boundaries  of  auroral  electron  precipi¬ 


tation  as  measured  by  DMSP/F2  satellite  with  the  inner  edge  of  the  plasmasheet 
in  a  very  large  statistical  sample,  we  find: 


1 


1)  Kp  orders  the  data  well,  particularly  in  the  evening  sector.  In 

a  given  MLT  sector,  the  latitude  of  the  equatorward  edge  of  the  precipitation 
decreases  linearly  with  Kp,  indicating  that  the  inner  edge  of  the  plasmasheet 
moves  earthward. 

2)  The  inner  edge  of  the  plasmasheet  is  earthward  of  the  Mauk-Mcllwain 
injection  boundary  for  Kp  >  2. 

3)  The  inner  edge  of  the  plasmasheet  is  not  symmetrical  about  the 
dawn-dusk  axis,  but  is  symmetrical  about  an  axis  whose  angle  to  the  dawn-dusk 
axis  depends  on  Kp,  reaching  a  constant  value  at  Kp  =  3.  That  is,  the  cross 
tail  electric  field  is  not  always  dawn  to  dusk  but  has  an  anti-sunward  component 
for  low  Kp  and  a  sunward  component  for  Kp  >  2. 

A)  The  inner  edge  of  the  plasmasheet  can  be  modelled  as  the  zero-energy 
Alfvdn  layer  in  a  magnetosphere  whose  large-scale  electric  field  is  obtained 
from  a  Volland-Stern  potential  with  y  >  2.  For  >  =  2,  the  Kp  dependence  of 
the  potential  is  found  to  be: 


.IB  R  ‘ 
o  E 


(3.9 


■OKp) 


10' 


L*-  sin  ( i  -4  )  -  - 
0  L 


Here,  4  is  dependent  on  1C  as  well,  but  is  not  easily  described  analytically. 
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1.  The  integral  flux  in  (cm  ster  sec)  (bottom  panel),  energy  flux 
in  keV  (cm^ster  sec)“l  and  average  energy  in  keV  (top  panel)  of 
precipitating  electrons  measured  by  the  DMSP/F2  satellite  passing 
over  the  south  pole  on  July  2,  1977.  These  are  plotted  as  functions 
of  universal  times  (in  seconds);  geographic  latitude  and  longitude 
of  the  sub-satellite  position;  corrected  geomagnetic  latitude  and 
longitude  of  the  satellite  position  projected  along  a  Jensen-Cain 
magnetic  field  to  100  km;  and  magnetic  local  time.  The  vertical 
lines  indicate  equatorward  and  poleward  boundaries  of  precipitating 
auroral  electrons. 


2.  Counts  in  selected  channels  of  the  SSJ/3  detector  plotted  as  a 
function  of  increasing  geographic  latitude  for  a  pass  over  the 
evening  side  auroral  zone  boundary. 

3.  As  in  Figure  2,  for  a  pass  over  the  morning  side  auroral  zone 
boundary . 

4.  Precipitating  electron  data  for  a  south  pole  pass  on  February  28, 
1978.  The  format  is  the  same  as  in  Figure  1. 

5.  Precipitating  electron  data  for  a  north  pole  pass  on  January  21, 

1977.  The  format  is  the  same  as  in  Figure  1. 

6.  Mean  values  and  standard  deviations  of  equatorward  boundaries  in 
each  line  plotted  as  a  functior  of  K  for  north  and  south 
poles  in  the  19-20  MLT  sector.  The  solid  line  results  from  a 
linear  regression  performed  using  individual  boundary  determinations. 
Linear  regression  results  are  summarized  for  given  Kp  ranges. 

7.  As  in  Figure  6  for  the  07-08  MLT  sector. 

8.  Equatorward  boundary  of  auroral  oval  as  a  function  of  MLT  for  high 
and  low  levels  of  magnetic  activity.  Data  marked  by  symbols 

0,  X,  A,  and  □  come  from  north  and  south  hemisphere  measurements 
of  DMSP/F2  (K  =  0,  5)  measurements  of  ISIS  (Bz  =  +  5>)  and 
DMSP  optical  Imagery  (K  =  0,  5)  respectively.  The  quiet  time 
oval  has  been  fit  to  an  offset  circle  (dashed  line). 

9.  Projections  of  the  circular  fit  of  Kp  =  0  boundaries  to  the 
ecliptic  plane  using  a  dipolar  projection  (dashed  line);  the 
Mead-Fairf ield  magnetic  field  model  (solid  line).  The  dash-dot 
curve  is  the  projection  by  means  of  a  Mead-Fairf ield  magnetic  field 
model  of  a  smooth  fit  to  the  individual  boundary  values  for  Kp  =  0. 

10.  A  comparison  of  the  Carpenter-Park  values  of  the  plasmapause 

at  dawn  to  the  inner  edge  of  the  plasmasheet  at  dawn  as  determined 
by  the  DMSP/F2  boundaries. 


Figure  11.  A  comparison  of  the  inner  edge  of  the  plasmasheet  as  determined 

by  the  DMSP/F2  boundaries  to  the  Mauk-Mcllvain  injection  boundaries 
for  various  Kp. 

Figure  12.  Volland-Stern  injection  boundaries  for  y  “  1  and  y  ■  2  rotated 
to  fit  the  inner  edge  of  the  plasmasheet  as  determined  by  the 
DMSP/F2  auroral  boundaries,  for  various  1C.  The  Mauk-Mcllwain 
injection  boundary  is  also  shown  (dashed  line). 
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COMMENTS  ON  "DIURNAL  VARIATION  OF  THE  AURORAL  OVAL  SIZE"  BY 


CHING  -I.  MENG 

by 

M.  S.  Gussenhoven, 1  D.  A.  Hardy,2  and  W.  J.  Burke1 

The  position  of  the  equatorward  boundary  (A  )  of  the  quiet  time  auroral 
oval  has  been  studied  by  Meng  (1979)  using  measurements  from  an  electron 
spectrometer  on  board  the  DMSP  (Defense  Meteorological  Satellite  Program) 
satellite  F2.  A  diurnal  vairation  of  ~4°  was  found  in  the  corrected  geomag¬ 
netic  latitudes  (CGL)  at  which  the  satellite  encountered  this  boundary.  The 
variation  was  attributed  to  a  physically  unexplained  30%  increase  in  the 
size  of  the  polar  cap  as  the  north  magnetic  pole  rotated  from  local  midnight 
to  local  noon.  Here  we  suggest  that  most  of  the  observed  variations  in  auroral 
boundaries  simply  reflect  diurnal  shifts  in  the  satellite's  trajectory  in  the 
magnetic  local  time  -  magnetic  latitude  (MLT-CGL)  frame  of  reference. 

The  DMSP  satellite  is  in  a  nearly  sun-synchronous  circular  orbit,  with  an 
inclination  of  98.8°,  centered  on  the  0700-1900  local  time  meridian.  Due  to 
the  offset  between  the  earth's  spin  axis  and  the  axis  of  the  magnetic  dipole 
the  orbit  has  significant  diurnal  variations  in  magnetic  frames  of  reference. 
This  is  illustrated  in  Figures  1  and  2  where  we  have  plotted  the  northern  and 
southern  high  latitude  portions  of  the  satellite's  MLT-CGL  trajectories  for 
21  January  1978  (the  quietest  day  in  Meng's  study).  These  trajectories  are 

1.  Physics  Department,  Boston  College,  Chestnut  Hill,  MA  02167 

2.  Space  Physics  Division,  Air  Force  Geophysics  Laboratory,  Bedford,  MA  01731 
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representative  of  those  followed  by  DMSP  during  the  1977-1978  winter  months. 

At  the  beginning  of  the  UT  day  the  north  polar  orbit  is  roughly  dawn  to  dusk 
crossing  the  noon-midnight  meridian  at  83°  CGL  in  the  night  hemisphere.  The 
most  poleward  excursions  of  the  orbit  occur  between  07  and  08  UT;  the  most 
equatorward  nighttime  excursions  are  found  near  1700  UT.  The  orbits  in  the 
southern  hemisphere  have  the  same  type  of  motion  in  the  opposite  direction. 

They  move  from  the  region  of  the  magnetic  pole  toward  the  cusp  (or  noon 
sector)  often  falling  outside  the  oval  entirely,  and  then  return  poleward. 

Since  the  oval  lies  at  higher  latitudes  in  the  noon  sector  than  in  the  midnight 
sector  there  are  fewer  possible  boundary  determinations  to  be  made  in  the  South 
Pole  data.  On  the  whole,  the  data  after  12  UT  on  a  given  day  is  unusable  for 
this  purpose  and  this  is  reflected  in  the  boundary  data  presented  by  Meng. 

Figure  3  illustrates  how  diurnal  variations  of  the  DMSP  orbit  give  rise 

to  systematic  variations  in  the  CGL  of  the  equatorward  auroral  boundary.  Here 

we  have  plotted  the  northern  hemisphere  values  of  A  detected  on  21  January  1978. 

h 

Meng  et  al.  (1977)  have  shown  that  the  poleward  boundaries  of  the  oval  (Ap) 
during  quiet  times,  as  marked  by  discrete  auroral  arcs,  are  well  represented 
by  offset  circles.  During  the  period  of  Meng's  study  emissions  from  the  oval 
were  below  the  detectability  levels  of  the  DMSP  optical  scanner.  We  have 
assumed  that  during  this  period  the  oval  is  latitudinally  narrow  and  that  the 
positions  of  A_  can  also  be  represented  by  offset  circles.  The  circle  in 
Figure  3  was  chosen  to  fit  the  morning  and  evening  values  of  Ap  for  the  first 
DMSP  pass  of  the  UT  day.  It  has  a  radius  of  18°  and  is  centered  at  86°  CGL  and 
0100  MLT.  The  circle  best  fitting  the  entire  data  set,  in  the  least  squares 
sense,  has  a  radius  of  18.5°  and  is  centered  at  86°  CGL  and  0115  MLT.  The 
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corresponding  southern  hemisphere  circle  is  of  19°  radius,  centered  at  -87°  CGL 
and  0120  MLT.  These  values  are  consistent  with  quiet  time  averages  reported 
by  Meng  et  al.  (1977).  As  seen  in  Figure  3  the  evening  oval  boundary  is  very 
stable,  lying  along  the  chosen  circle.  The  morning  oval  boundary  is  close  to 
the  plotted  circle  but  shows  much  greater  variability.  This  is  consistent 
with  results  from  a  large  set  of  DMSP  measurements  reported  by  Gussenhoven  et  al. 
(1979)  who  sv  w  that:  (A)  The  morning  equatorward  boundary  is  much  more 
variable  than  the  evening  boundary  for  a  given  magnetic  activity.  (B)  For 
extended  periods  of  quiet  magnetic  activity  the  oval  does  not  become  entirely 
stable,  but  can  fluctuate  greatly  in  its  poleward  boundaries  and  to  a  lesser 
extent  in  its  equatorward  boundaries. 

Predicted  UT  modulations  of  A„  found  by  combining  offset  circles  and  DMSP 
orbital  variations  for  the  northern  and  southern  hemispheres  are  given  in 
Figures  4  and  5  respectively.  The  upper  (lower)  traces  represent  the  afternoon¬ 
evening  (morning)  equatorial  boundaries.  Observed  values  of  A  ,  as  determined 
by  Meng,  are  plotted  as  triangles.  Both  the  afternoon  and  morning  southern 
hemisphere  as  well  as  the  evening,  northern  hemisphere  observations  are  in 
remarkable  agreement  with  their  predicted  estimates.  The  morning,  northern 
hemisphere  measurements  show  less  agreement.  The  overall  differences  would  be 
less  had  we  used  the  least  squares  circle  as  our  predictor.  It  should  also 
be  noted  that  morning  and  evening  values  tracked  each  other  better  on  the  other 
days  presented  by  Meng  than  they  did  on  the  first  half  of  21  January.  Thus, 

on  the  other  days  the  predicted  and  observed  values  of  A  would  be  in  better 

b 

agreement  than  in  the  case  presented  here. 
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We  therefore  conclude  that  a  very  large  portion  if  not  all  of  the  diurnal 
auroral  boundary  variation  reported  by  Meng  can  be  explained  by  variations 
in  the  satellite  orbit  across  an  oval  represented  by  an  offset  circle.  We 
suggest  that,  in  order  to  determine  the  universal  time  variation,  either  a 
correction  for  the  effect  of  orbit  change  in  local  time  could  be  made,  or  that 
a  much  greater  restriction  in  magnetic  local  time  be  imposed. 
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Figure  1.  The  orbital  paths  of  the  DMSP-F2  satellite  over  the  North  pole 


Figure  2 
Figure  3 


Figure  4 


for  January  21,  1978.  The  orbits  are  plotted  in  corrected 
geomagnetic  latitudes  and  magnetic  local  time.  The  approximate 
universal  time  for  each  polar  crossing  is  given.  The  quiet  oval 
is  represented  by  an  offset  circle. 

Same  as  Figure  1,  but  for  South  pole  crossings. 

The  boundaries  for  the  northern  auroral  zone  as  determined  by  Meng, 
plotted  in  corrected  geomagnetic  latitude  and  magnetic  local  time 
for  January  21,  1978,  along  with  the  offset  circle  used  to  predict 
the  diurnal  variation.  The  small  cross  marks  the  geomagnetic  pole, 
the  large  cross  the  center  for  the  offset  circle. 

The  solid  lines  show  the  expected  dirunal  variation  in  the  equator- 
ward  boundary  of  the  auroral  zone  in  the  northern  hemisphere  for 
both  morning  and  evening  sectors  when  the  variation  in  the  orbit 
has  been  taken  into  account  and  when  the  oval  has  been  represented 
by  an  offset  circle.  The  evening  boundaries  are  shown  in  the  upper 
panel  (for  latitude  marks  on  the  left)  and  the  morning  boundaries 
in  the  lower  panel  (for  latitude  marks  on  the  right).  The  boundaries 
as  determined  by  Meng  are  plotted  as  closed  (evening  side)  and  open 
(morning  side)  triangles. 


Figure  5.  Same  as  for  Figure  3,  but  for  the  South  pole. 
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